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Q ■ ABSTRACT 

<N ' 

Aims. We study the resolved stellar populations and derive the star formation history of the SDIG, a gas-rich dwarf galaxy member 
■ of the NGC7793 subgroup in the Sculptor group of galaxies. 

Methods. We construct a color-magnitude diagram using archival observations from the HST/ACS in order to examine the stellar 
content of SDIG, as well as the spatial distribution of stars selected within different stellar evolutionary phases. We derive the star 
formation history of SDIG using a maximum-likelihood fit to the color-magnitude diagram. 

Results. The color-magnitude diagram shows that SDIG contains stars from 10 Myr to several Gyr old, as revealed from the main 
sequence, blue loop, luminous asymptotic giant branch, and red giant branch stars. The young stars with ages less than ~250 Myr 
show a stellar spatial distribution confined to the central regions of SDIG, and additionally the young main sequence stars exhibit an 
off-center density peak. The intermediate-age and older stars as traced by the red giant branch stars are more spatially extended. SDIG 
is dominated by intermediate-age stars with an average age of 6.4+j Gyr. The average metallicity inferred from the CMD modelling 
is [M/H] as -1.5 dex. SDIG has a star formation history consistent with a constant star formation rate, except for ages younger than 
s;200 Myr. The lifetime average star formation rate is 1.3^Sj xl0~ 3 Mq yr -1 . More recently than 100 Myr, there has been a burst of 
star formation at a rate ~2-3 times higher than the average star formation rate. The inferred recent star formation rate from CMD 
modelling, 2.7(±0.5) xl0~ 3 Mq yr -1 , is higher than inferred from the Ho- flux of the galaxy; we interpret this to mean that the upper 
end of the initial mass function is not being fully sampled due to the low star formation rate. Additionally, an observed lack of bright 
blue stars in the CMD could indicate a downturn in star formation rate on 10 7 -yr timescales. A previous star formation enhancement 
appears to have occurred between 600-1100 Myr ago, with amplitude similar to the most recent 100 Myr. Older bursts of similar 
peak star formation rate and duration would not be resolvable with these data. The observed enhancements in star formation suggest 
that SDIG is able to sustain a complex star formation history without the effect of gravitational interactions with its nearest massive 
galaxy. Integrating the star formation rate over the entire history of SDIG yields a total stellar mass equal to 1J7^J1 xlO 7 Mq, and 
a current V-band stellar mass-to-light ratio equal to 3.2 Mq / Lq. 

. 
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t—( ' 1. Introduction of individual SFHs is valuable to uncover and understand the 

£N| . intrinsic and environmental differences affecting their evolution 

, Resolved stellar populations and star formation histories (SFHs) ( e g \jy e i sz e ( a \ |2Q1 lbt 
• • ■ provide a powerful tool to understand galaxy formation and evo- 

lution in all scales (e.g., Cignoni & Tosi 2010). Detailed SFHs ° ne such interesting galaxy to focus on is the Sculptor 
of Local Group (LG) dwarf galaxies have shown that these are Dwarf Irregular Galaxy (SDIG). First discovered by Laustsen 
complex systems (e.g., MateoEM! Tolstoy, Hill & Tosi 2009). et aL <EZ3' SDIG is a low-luminosity (Heisler et al. [T993, 
£ 1 The complexity of SFHs observed in LG dwarfs has been ex- f s - nch £^L f T 11Te g ular ^ ^vsky Falgarone & 

plored by means of numerical simulations and the results show Lequeux[T977l Lequeux & West 1981 ; Cote et aLHSSZ). Atomic 
that their diversity can be understood by invoking the action hydrogen observations show that SDIG has a smooth distribution 
of ram pressure stripping, tidal stirring, and/or cosmic ultravi- wlthout the presence of structures such as shells, holes or bub- 
olet (UV) background radiation upon gas-rich dwarf galaxies bles ^ Cote et d-EED- SDIG contains no H II regions and this is 
(e.g., Mayer [2010, Kazantzidis et al. 201 1 , Nichols & Bland- suggestive of a quiescentphase with very lowcurrent star forma- 
Hawthom[20TTl and references therein). From an observational tlon rate < SFR; Miller[T996l Heisler etaL[T997l Skillman, Cote 
point of view, the comparison of SFHs between dwarf galaxies & Mlller ™D. Bouchard et al. 420® detect two Hor point 
in the LG and outside the LG shows that the former are represen- sou f ces and estlmate a ver y low current (3.8xl0 5 M 
tative of their kind (Weisz et al. l2TJTTa1 l. so that a deeper study y r > ; however ' Lee et aL ™) detect far-ultraviolet (FUV) and 

near-ultraviolet (NUV) emission that appears to have a clumpy 

morphology and suggest main-sequence B stars as the source of 

* Based on observations made with the NASA/ESA Hubble Space ■ ■ 

_ , . . , . , , , . , _ ™ i this emission. 

Telescope, obtained trom the data archive at the Space Telescope 

Institute. STScI is operated by the association of Universities for With a distance of 3.2 Mpc and based on its radial veloc- 

Research in Astronomy, Inc. under the NASA contract NAS 5-26555. ity, SDIG is a member of a loose galaxy triplet (Karachentsev 
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et al. 120061 1. of which the other two galaxy members are 
NGC7793 at a distan ce of 3.61+0.53 Mpc (Vlajic, Bland- 
Hawthorn & Free man 120111) and UGCA442 at a distance of 
3.8+° * Mpc (Mould 2005J. This triplet, the NGC 7793 subgroup 
(Karachentsev et al. 120031 ) lies at the far side of the Sculptor 
group. SDIG therefore lies near the limit of the ACS Nearby 
Galaxy Survey Treasury galaxy sample (ANGST; Dalcanton et 
al. [2009), so a study of the SFH of this subgroup represents 
an opportunity to learn more about the Sculptor group, a low- 
density and extended filamentary structure with several embed- 
ded subgroups (e.g., Jerjen, Freeman & Binggeli 1998 ). The fil- 
amentary structure of the Sculptor group extends from 1 .5 Mpc 
to 4 Mpc, and points to an unvirialised state, with galaxies still 
falling in (Jerjen, Freeman & Binggeli 1998 , 2000 1 Karachentsev 
2005 ). The early evolutionary state and the low density of the 
Sculptor group provide an ideal environment in which its galaxy 
members may be considered as still evolving in isolation. 

In the present work, we examine the resolved stellar popula- 
tions and model the color-magnitude diagram (CMD) of SDIG 
using archival Hubble Space Telescope (HST) observations. The 
structure of this work is as follows. In §2 we present the obser- 
vations and photometry, including photometric error analysis. In 
§3 we show the CMD and discuss the stellar content of SDIG, in- 
cluding an examination of the spatial distribution of stellar pop- 
ulations selected in several stellar evolutionary phases. In §4 we 
briefly review the CMD modelling technique. In §5 we show 
the results for the SFH of SDIG. We summarise and discuss our 
findings in §6. 

2. Observations and photometry 

We use archival observations taken with the Advanced Camera 
for Surveys (SNAP 9771; Karachentsev et al. 120061 ) with a 
total exposure time of 900 sec in the F814W-band, and 
1200 sec in the F606W-band. We retrieve the pre-reduced im- 
ages through the ST-ECF Hubble Science Archive and per- 
form stellar point source photometry using the ACS module 
of DOLPHOT, a modified version of the photometry pack- 
age HSTphot (Dolphin 2000). We apply point source pho- 
tometry simultaneously on all the individual, calibrated and 
flat-fielded FLT images, and we use as a reference image 
the drizzled image in the F814W filter. We proceed with the 
photometry reduction steps as described in the ACS module 
of DOLPHOT. DOLPHOT provides magnitudes in both the 
ACS/WFC and the Landolt UBVRI photometric systems and 
we choose to use the ACS / WFC filter system for our study. The 
final photometric catalogue includes 1 1363 stars selected to have 
S/N > 5, sharpness |sharpF606w+sharpp8i4wl <L and crowding 
(crowdF606w+crowdp8i4w) <1 (e-g-, Williams et al. 120091 ). 

We conduct artificial star tests in order to quantify the pho- 
tometric errors and completeness of our data, using the utilities 
provided within DOLPHOT. To that end, we run artificial star 
tests with 10 5 stars per ACS /WFC chip, with range in F606W 
from 21 mag to 30 mag, and in (F606W-F814W) from -1.5 mag 
to 3.5 mag. In order to better quantify the photometric errors 
in the faintest magnitudes, which we need for the CMD mod- 
elling, we run complementary artificial star tests with an addi- 
tional number of 2xl0 5 stars per ACS/WFC chip and cover- 
ing magnitudes fainter than 25 mag. The high number of stars 
inserted during the artificial star tests is not an issue for self- 
induced crowding, since each star is inserted and photometered 
one at a time (Dolphin 2000, see also Perina et al. 2009 ). 

The photometric errors derived using the artificial star tests 
are shown in Fig. Q] The mean photometric error and standard 
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Fig. 1. Photometric errors derived from artificial star tests and es- 
timated as the difference between output and input magnitudes, 
as a function of the output magnitude in the F814W filter (upper 
panel) and in the F606W filter (lower panel). 
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Fig. 2. Completeness factors as a function of the F606W-band 
magnitude, upper panel, and F814W-band magnitude, lower 
panel. The red arrow in each panel indicates the 50% complete- 
ness factor. 



deviation is 0.03+0.1 mag in F814W, and 0.02+0.09 mag in 
F606W. In Fig. [2] we show the completeness factors as a func- 
tion of magnitude. The 50% completeness factor occurs at the 
magnitude F606W =~26.9 mag, and F814W =~25.9 mag. 

Table[T]lists the basic properties for SDIG. Each row shows: 
(1) the galaxy type; (2) and (3) equatorial coordinates of the 
field center (J2000.0); (4) the reddening adopted from Schlegel 
et al. ( 119981 ); (5) the galactic foreground extinction in the F606W 
and F814W filters, computed using the extinction ratios for a 
G2 star listed in Sirianni et al. (2005, their Table 14); (6) the 
magnitude of the tip of the red giant branch (TRGB) we derive 
from the data; (7) the distance modulus we derive from the data; 
(8) the derived distance in Mpc; (9) the absolute V-band mag- 
nitude adopted from Georgiev et al. ( 120091 ), who use the dis- 
tance of Karachentsev et al. 2006; (10) the H i mass adopted from 
Koribalski et al. (2004); (11) and (12) the magnitude-weighted 
coordinates of the center, xc and yc, respectively, which we 
compute from the data (see §3); (13) and (14) the ellipticity and 
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Table 1. Properties of SDIG. 



Quantity 


Value 


Type 


dlrr 


RA (J2000.0) 


00*08'" 13.36 s 


Dec (J2000.0) 


-34°34'42.0" 


E(B - V) (mag) 


0.012 


^Feoew^Fsww (mag) 


0.034, 0.022 


F814W T rgb (mag) 


23.44+0.05 


(m - M) (mag) 


27.51 ±0.06 


Distance (Mpc) 


3.2+0.1 


M v (mag) 


-11.87 


Mm (M©) 


2.5xl0 7 


x c (arcsec) 


105.1±1.7 


yc (arcsec) 


114.5+1.3 


e 


0.05 


PA (degrees) 





r cff (arcsec) 


29.4 ± 1.6 



22 - 




-0.5 0.5 1 1.5 2 
F606W - F814W 



Fig. 3. Color-magnitude diagram of SDIG. The gray dashed line 
shows the 50% completeness factors. The error-bars correspond 
to the photometric errors derived using artificial star tests. 

position angle adopted from Kirby et al. (120081) : (15) the effec- 
tive radius adopted from Kirby et al. J20081 I. 

3. Color-magnitude diagram and stellar content 

The CMD of SDIG shown in Fig.[3]reveals the presence of sev- 
eral stellar evolutionary features. The most prominent feature is 
the red giant branch (RGB). The RGB indicates that stars with 
ages from ~1.5 Gyr and older are present (e.g, Salaris, Cassisi 
& Weiss 2002). Intermediate-age stars, with ages between 1 Gyr 
and 10 Gyr, are revealed in the presence of luminous asymptotic 
giant branch (AGB) stars, while the young main-sequence (MS) 
and blue-loop (BL) stars indicate the presence of stars with ages 
younger than 1 Gyr. Therefore, SDIG contains stars spanning a 
wide range in age, from ~Myr to older ages. 

In the following analyses and unless otherwise noted, we 
use the TRGB magnitude and distance modulus we measure 
in the ACS / WFC filter system. We derive the distance modu- 
lus for SDIG using the absolute F814W-band magnitude of the 
TRGB computed through the calibration of Rizzi et al. (120071) : 



M F814W = -4.06 + 0.20 x [(F606W - F814W) - 1.23], and the 
F814W-band magnitude of the TRGB computed through the 
Sobel-filtering technique on the F814W-band luminosity func- 
tion (e.g., Lee, Freedman & Madore 1993; Sakai et al. 1996). 
We find a TRGB F814W-band magnitude of 23.44+0.05 mag, 
while measuring a dereddened color at the TRGB magnitude 
level of 1.05±0.03 mag, we derive M F8 i 4w = -4.09+0.02 mag. 
The quoted uncertainties take into account the photometric er- 
rors from the artificial star tests and the width of the Sobel filter 
response. These values along with the extinction listed in TableQ] 
give a distance modulus of 27.51 ±0.06 mag, which translates to 
a distance of 3.2+0.1 Mpc. The quoted uncertainties take into 
account the width of the Sobel filter response for the TRGB 
magnitude detection, uncertainties related to the TRGB calibra- 
tion, and the photometric errors. Our derived distance modulus 
and distance are consistent with the values of 27.53 mag and 
3.21 Mpc, respectively, found in Karachentsev et al. ((2006), who 
use the same method on the transformed V and I filters and as- 
sign an uncertainty to their distance of 8%. We list the derived 
F814W-band magnitude of the TRGB, distance modulus, and 
distance in Table [TJ 

Fig.[4]shows the CMD of stars distributed within three ellip- 
tical radius annuli: within 1 effective radius, r e ft-; between 1 r e e 
and 2 r e s\ and beyond 2 r e g\ The elliptical radius is defined 
as: r — [x 2 + y 2 /(l - e) 2 ] 1 ^ 2 , where x and y are the distances 
along the major and minor axis, respectively, and e is the el- 
lipticity. We adopt the ellipticity, position angle, and effective 
radius from Kirby et al. (2008), listed in Table [TJ Also listed in 
TableQJare the xc and yc coordinates of the magnitude-weighted 
center of SDIG, which we compute from our data using the 
equation: xc = £(mag, x x,)/£mag ; , where mag, corresponds 
to the F814W-band magnitude of each star, and x, to their x (or 
y for the yc) coordinates, in pixels. We use stars with magni- 
tudes brighter than the magnitude with 50% completeness fac- 
tors. In the right panel of Fig. [4] we show with the open circles 
the Galactic foreground contamination, which was estimated us- 
ing the TRILEGAL cod e (Gir ardi et al. [20051 Vanhollebeke, 
Groenewegen & Girardi 120091 ). We expect a total of 58 fore- 
ground stars in the same color and magnitude location of the 
SDIG stars, which translates to less than 0.5% of the total num- 
ber of SDIG stars. In the left panel of Fig.[4]we overplot Padova 
isochrones (Marigo et al. 120081 Girardi et al. [2010), with ages 
between 20 Myr and 150 Myr. Fig. [4] indicates that the bulk 
of the young MS stars are confined within the central parts of 
SDIG. 

We examine the cumulative distribution functions of the 
MS, BL, luminous AGB, and RGB stars, versus the elliptical 
radius, r, defined in the preceding paragraph. We select those 
MS stars with their color ranging from -0.5 mag to -0.1 mag 
and their F814W-band magnitude ranging from 23.8 mag to 
25.25 mag; the BL stars are selected with their color ranging 
from -0.05 mag to 0.35 mag and with their F814W-band mag- 
nitude ranging from 21 mag to 25 mag. The RGB stars are se- 
lected between the magnitude of the TRGB and 1 mag below, 
while the color selection ranges such as to follow the slope of 
the RGB. The luminous AGB stars are selected to be brighter 
by 0.15 mag than the TRGB ( Armandroff et al. [19931 and to lie 
within 1 mag above (Itrgb - 0.15) mag, and within the color 
range of 1 < (V - I)o < 3.5 (mag). We count a number of 80 lu- 
minous AGB stars, which translates into a fraction, i.e. number 
of luminous AGB stars versus the number of RGB stars within 
1 mag below the TRGB magnitude, N(AGB) / N(RGB)=13%, 
comparable to that of Sculptor group dwarf galaxies (Lianou et 
al. I20121 ). The cumulative distribution functions in Fig. [5] show 
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Fig. 4. Color-magnitude diagrams within three elliptical radius annuli. The first annulus corresponds to an elliptical radius equal to 
1 r e ff , the second between 1 and 2 r e ff, and the third beyond 2 r e g, shown from left to right, respectively. In the left panel, we overplot 
Padova isochrones with ages of 20 Myr, 40 Myr, 60 Myr, 80 Myr, and 150 Myr, and for all isochrones a Z metallicity of 0.001. 
Also shown in the left panel with the error bars the photometric errors. In the right panel, we show with open circles the galactic 
foreground contamination. 




RGB: thick solid 
AGB: thick dashed 
BL : thin solid 
MS: thin dashed 
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Fig. 5. Cumulative distribution functions for the selected stellar 
populations. The vertical line shows the effective radius, r c yy. 



that the MS and BL stars are confined to the central parts of 
SDIG, while the luminous AGB and the RGB stars are more spa- 
tially extended. Performing a two-sided Kolmogorov-Smirnov 
(K-S) test between each of the cumulative distributions for the 
MS, BL, and luminous AGB stars with the cumulative distribu- 
tion of the RGB stars, we derive results consistent with spatially 
separated populations at the 99% confidence level in all cases. 

Comparing the spatial distribution of the luminous AGB 
stars with that of the RGB stars, in the stellar density maps 
shown in Fig. [6] the former are more confined in the inner parts 
of SDIG. The selected BL stars are confined to the central parts 
of the dwarf when we compare their spatial distribution with the 
one of the luminous AGB and RGB stars. The stellar density map 
of the young MS stars shows a population concentrated towards 
the central parts of SDIG, while in the same time its spatial lo- 
cation is not coincident with the very center of SDIG, having 
an offset of 8.2 arcsecs, or 130 pc at a distance of 3.2 Mpc. We 
note that a similar characteristic of off-center local sites of young 
star formation characterise the majority of the transition-type 
dwarfs in the Sculptor group (Jerjen & Rejkuba 2001 ; Lianou 
et al. 120121 1. while the same is also observed for the young stel- 
lar populations of other LG dwarf galaxies (e.g, Dohm-Palmer 
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Fig. 6. Gaussian-smoothed stellar density maps for stars selected 
in the MS, BL, luminous AGB, and RGB phases, expressed in 
terms of angular offsets from the dwarf's center. The density 
maps are color coded such that black corresponds to a density of 
0.16 stars per sq. arcsec and the density declines toward lighter 
grey. The ellipse has a major axis equal to 1 r e g , while the crosses 
emphasise the location of the center of the dwarf. The red arrows 
in the panel of the RGB population indicates the directions of 
North and East, and is the same for all panels. 



et al. [T997l Gallagh er eta l. [T998l Mateo fl998l Cole et al-I1999; 
McConnachie et al. 120061 1. 



4. Color-magnitude diagram modelling 

The basic idea using the modelling of the CMD to uncover the 
evolutionary history of a galaxy is that its composite stellar pop- 
ulations can be regarded as a linear combination of simple stel- 
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lar populations of different ages and metallicities. Thus, the SFH 
can be obtained by comparing the observed features in the CMD 
with those in a model CMD. There are several features in a CMD 
that can be used as age and metallicity indicators (e.g., Gallart, 
Zoccali & Aparicio 2005 ). Synthetic populations are constructed 
using a set of stellar evolution models, adopting an initial mass 
function (IMF) and applying a star formation and a chemical en- 
richment law (e.g., Gallart, Zoccali & Aparicio l2005l Cignoni & 
Tosi[2pi0). In practice, the density of stars within a defined color 
and magnitude grid along the observed CMD is compared with 
the density of the synthetic stars drawn from the model CMD 
(for a review see Tolstoy, Hill & Tosi 120091 Aparicio|2002), and 
the best model is adopted as the SFH. 

The approach taken here was developed by Cole to study the 
stellar populations of dlrr galaxies in the LG (see, e.g., Skillman 
et al. l2003bl Cole et al. 120071 Cignoni et al. 120121 . As with any 
such approach, the power of the method is greatly reduced by 
lack of photometric depth (e.g., Weisz et al. 201 lb ). Detailed in- 
ferences about the SFH will be restricted to age ranges for which 
the age-sensitive MS, BL, and to a lesser extent AGB stars are 
visible. It is still possible to place constraints on long-term av- 
erage SFRs and mean metallicities based on inferences from the 
RGB and the regions of the CMD that are badly affected by in- 
completeness. However, these inferences are by necessity more 
model-dependent and of lower precision, and the resulting un- 
certainties are dominated by systematic effects. 

The SFH-fitting starts with theoretical isochrones interpo- 
lated to a fine grid of age and metallicity in order to create a 
synthetic CMD with no gaps. The most recent models from the 
Padova models are used. The synthetic CMDs are binned in age 
and metallicity (Z) to avoid "overfitting" noise in the CMDs. 
The bins are initially spaced by a default value of Alog(age) = 
0.10, but adjacent bins are merged when lack of information in 
the CMDs leads to instability in the solutions. In this case, ages 
greater than 1 Gyr are only represented on the RGB and AGB, 
and the age resolution is decreased to Alog(age) = 0.30 (factor 
of two). There are not many stars younger than 30 Myr in the 
CMD, so the youngest age bin spans the entire range from 4- 
30 Myr as a way to populate the CMD bins sufficiently to draw 
SFH inferences. 

The CMD is divided into a regular grid of color-magnitude 
cells, and the expectation value of the number of stars in each 
cell for a SFR of 1 M Q yr _1 is calculated from the isochrones. 
Several parameters are assumed to be fixed during the solution: 
the distance modulus and reddening, IMF, fraction of binaries, 
and the binary mass ratio distribution function. The adopted IMF 
is from Chabrier (2003), and the binary fraction and mass ra- 
tios are parameterized based on Duquennoy & Mayor (1991) 
and Mazeh et al. ([1992). We take 35% of stars to be single and 
the rest to be binary. The binaries are divided into "wide" and 
"close" binaries in a 3:1 ratio; the secondaries in the wide sys- 
tems are drawn from the same IMF as the primaries, but in the 
close systems the secondary masses are drawn from a flat IMF. 
The distance and reddening are initially constrained to the values 
given in Table 1, but are varied if the resulting synthetic CMDs 
are mismatched to the data. 

No age-metallicity relation is explicitly assumed, but a range 
of metallicities at each age is allowed, constrained by the color 
range of the data. Because of the shallow photometric depth of 
the data, there are few strong constraints on the metallicity, and 
those that exist are model-dependent. The RGB is affected by 
age-metallicity degeneracy, but indicates a range of metallicities 
between 0.0004 < Z < 0.002. These numbers are consistent with 
the Padova AGB star tracks, and the color separation between the 



BL stars and MS. There is no metallicity constraint from H n re- 
gions, but the inferred values are consistent with the metallicity- 
luminosity relationship for dlrr galaxies in Lee, Zucker & Grebel 
( 12007b (see discussion below). 

The IMF-weighted, color-magnitude binned isochrones are 
convolved with color and magnitude error distributions derived 
from the artificial star tests in order to create synthetic CMDs 
with the same properties as the data. Linear combinations of the 
synthetic CMDs are created by the fitting routine, which uses 
a maximum likelihood test based on the Poisson distribution of 
counts in each bin (Cash 119791 ) to find the combination that is 
most likely to produce the observed CMD. We use a simulated 
annealing technique to find the best fit while avoiding false, local 
maxima in likelihood space. The SFH errorbars are calculated by 
successive perturbation of each age-metallicity bin; because the 
total number of stars is fixed by the observations, a decrease in 
SFR in one bin results in an increase in adjacent bins; this gives 
anticorrelations between adjacent bins. In general, the code is 
driven by the most populous cells of the CMD to its preferred so- 
lution; because most stars are observed on the MS and low-mass 
stars vastly outnumber those of higher mass, the most significant 
cells in determining the SFH are frequently below the 50% com- 
pleteness threshold, and the solution depends sensitively on the 
artificial star test results. 

Our data are not deep enough to reach the oldest MS turn- 
offs, and due to the age-metallicity degeneracy on the RGB, the 
solution of the CMD modelling for ages higher than ~1.5 Gyr 
bears larger uncertainties. The observed and modelled CMDs for 
SDIG are shown in Fig. [7] From the CMD modelling, we derive 
a distance modulus of 27.54 mag, while the reddening required 
is E(B-V)=0.04. The higher reddening, as compared to the value 
from Schlegel et al. ( 1998 ), is consistent with the one that Heisler 
et al. d 19971) derive for the very central parts of SDIG. 

The best-fit solution under the restrictions of the fit proce- 
dure matches the mean stellar density and colors of the major 
stellar sequences well, but there are some notable differences be- 
tween the simulated CMD and the data. The MS and BL in the 
data are blurred, with no clear color gap between them; however, 
in the models these features are well-separated. Presumably this 
is due to differential reddening, with a tendency for young stars 
to be associated with dust clouds and circumsteilar material. 

5. Star formation history 

We show the best-fit SFH for SDIG in Fig.g] Overall, the best-fit 
SFH is characterised by the following features: 

- An average age for formation of stars of (r) — 6.4+ j'jj Gyr and 
a total astrated mass (time-integrated SFR) of 17.7x10 6 Mq. 

- The average metallicity is inferred to be Z «0.0006 ([M/H] 
w — 1.5 dex), with some tendency to increase with time, as 
the BL stars are more consistent with a metallicity closer to 
Z = 0.001. 

- A long-term SFH that is consistent with a constant SFR over 
several Gyr. The average SFR is 1.3+°3 xl0~ 3 MQyr 1 . Any 
bursts or fallow periods are of short enough duration and/or 
amplitude that they do not strongly perturb the mean SFR 
averaged over factor of two periods in age. This is broadly 
consistent with the typical dlrr galaxy found in Weisz et 
al. (1201 lbk however the HST SDIG observations analysed 
here were not deep enough to reveal factor of 2-3 variation 
in SFR over sub-Gyr intervals at intermediate or old ages, 
and these are not excluded by the CMD analysis. 
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Fig. 7. Left panel: logarithmically-scaled, binned CMD used in the maximum-likelihood fit for SFH; right panel: best-fit model. 
Without differential reddening, there is clear separation between the MS and BL stars, which are mingled together in the data. The 
AGB stands out as being poorly-fit; the mismatch occurs because of the large number of stars fainter than m814 x26, which are part 
of the same stellar population as the AGB stars. 



- The age resolution improves dramatically for ages less than 
~800-1000 Myr, but strong variations over intervals of a few 
hundred Myr are not observed, with two exceptions: 

1 . Over the past 100 Myr the mean SFR increases by a fac- 
tor of «2 over its long-term average value, but then de- 
clines slightly. Variations of this duration and amplitude 
would be completely hidden by photometric incomplete- 
ness and age-metallicity degeneracy for ages older than 
~1 Gyr; this behaviour may be typical of the history of 
low-luminosity dlrrs (e.g., Cole 2010). Because the frac- 
tion of evolved stars is low, the derived SFR is suscep- 
tible to uncertainties in the modelling of incompleteness 
and possible variations in the IMF. Stars in this age range 
would contribute to the SFR as measured from FUV data, 
but not to the Ha luminosity owing to the low absolute 
SFR with a value of -2.5+^ xl0~ 3 Mq yr~ x from 0- 
10 Myr. 

2. A similar-amplitude SFR enhancement is inferred to 
have taken place from 600-1 100 Myr ago. The younger 
stars in this interval are directly observable in the CMD 
as BL stars at F606W-band magnitude m 26.5-27 mag on 
the blue side of the RGB. The older stars are confused 
with the RGB and highly subject to incompleteness, re- 
sulting in the large uncertainties in the SFR at this age 
range. Because of the high incompleteness, the best-fit 
SFR has high uncertainty (reflected in the model over- 
production of AGB stars associated with this age range). 
Based on this data alone the balance of evidence favours 



the existence of a <500-Myr long episode of increased 
SFR ss 1 Gyr ago, but is consistent with a flat SFH given 
the small number of stars and the statistical nature of the 
problem. As above, this type of event would be unde- 
tectable at older ages based on this data set alone. 

Integrating the SFR over the entire history gives a total as- 
trated stellar mass M*= 1.77+° ^ xlO 7 Mq. For ages less than 
1 Gyr, we estimate a stellar mass M^ioyr = 1.68^ ® xlO 6 Mq, 
which consists a fraction of ~10% of the total stellar mass. 
Because the maximum mass of surviving stars decreases 
with age, only ^85% of the total astrated mass is still 
present in the form of stars, i.e., the current stellar mass 
M*, cun - ent = 1.50+Jj ^ xlO 7 M©. Kirby et al. (120081 derive a total 
stellar mass equal to 1.6xl0 7 Mq, using deep near-infrared ob- 
servations that trace the old stellar populations and assuming an 
H-band stellar mass-to-light ratio of 1 . The agreement between 
our current stellar mass and the total stellar mass of Kirby et 
al. (120081 1 is very good. 

The current stellar mass is comparable to the Hi mass 
in SDIG (see Table [Q Koribalski et al. 12004b with a ratio 
of Mm / M* = 0.6. We compute the baryonic gas fraction, 
Mga S / Mbaryonic, and compare it with the baryonic gas frac- 
tion of the ANGST sample galaxies (Weisz et al. 1201 lbl their 
Fig. 9). The baryonic gas fraction of SDIG is estimated equal to 
M gas / Mbaryonic =0.7, assuming a gas mass M g a S =1.4xMhi and 
a baryonic mass Mbaryonic =M g a S +M*. These values are listed in 
Table [2] Placing SDIG in Fig°9 of Weisz et al. ( 1201 lbl l locates 
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Fig. 8. Star formation rate as a function of age for SDIG based 
on the maximum-likelihood fit, shown in the upper panel with 
a linear scale on age, and in the lower panel with a logarithmic 
scale with age. Note in the upper panel the break in age scale at 
750 Myr that reflects the decreasing time resolution of the CMD 
with increasing age. 



it in the upper left region occupied by dlrrs, with SDIG slightly 
above ESO321-G014. 

Based on our current stellar mass estimate and the V-band 
absolute magnitude listed in Table [T] we compute a V-band stel- 
lar mass-to-light ratio equal to M* / Ly = 3.2 M© / Lq. This 
is higher than transition-type LG galaxies such as LGS3 and 
DDO 216 (Hunter, Elmegreen & Ludka [2010l l. Using the B-band 
magnitude listed in Karachentsev et al. (2004; column (6) in their 
Table 1), we compute a B-band stellar mass-to-light ratio equal 
to M* / L B =1 .4 M Q / Lq . The B-band stellar mass-to-light ratio 
of SDIG is lower than the average B-band stellar mass-to-light 
ratio (equal to 2.4 Mq / Lq) of the sample of star-burst dwarf 
galaxies listed in McQuinn et al. (2010J, but is higher than the 
average value of ~ 1.1 Mq / Lq for dlrrs with -11 mag < Nig < 

-13 mag (Weisz et al.EfllB- 

SDIG is a late-type dwarf galaxy and the dominant popula- 
tion has intermediate ages, with an average age of 6.4^]^ Gyr. 
Thus, the age-metallicity degeneracy on the RGB makes it diffi- 
cult to derive individual stellar metallicities. On the other hand, 
the mean stellar metallicity of mixed stellar populations re- 
mains rather robust against the age-metallicity degeneracy on 
the RGB (Lianou et al. 2011). Therefore, we use the average 
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Fig. 9. CMD focused on the RGB stars. The grey circles indicate 
the RGB stars used to compute the fiducial line, with the lat- 
ter shown with the open black circles. The solid red line shows 
the Dartmouth 6.4 Gyr, solar-scaled isochrone with metallicity 
-1.6 dex in [M/H] that best matches the RGB fiducial, while 
the black dotted line shows the Padova 6.4 Gyr isochrone with 
metallicity -1.5 dex in [M/H]. Towards fainter magnitudes, the 
two isochrones overlap. 



age of the stellar populations in SDIG in order to constrain its 
mean stellar metallicity. To that end, we use the fiducial of RGB 
stars, with magnitudes ranging from the TRGB to 1 mag fainter 
and shown in Fig. [9] with the open circles, to which we over- 
plot the best-matched isochrone from the Dartmouth (Dotter et 
al. 2007 , 2008) and Padova Stellar Evolution Databases, with an 
age equal to the average age of the stellar population in SDIG. 
Here, we use the colour excess and distance modulus we derive 
from the SFH modelling (see §5). The best-matched isochrones 
are shown in Fig. [9] with the red solid line and black dotted 
line, respectively, for the Dartmouth and Padova isochrones. 
The metallicity of the best-matched Dartmouth isochrone is 
— 1.6 dex, while the metallicity of the best-matched Padova 
isochrone is —1.5 dex. The metallicity of the Padova isochrone 
that best matches the RGB fiducial is consistent with the aver- 
age metallicity we derive from the SFH-modelling. The inferred 
stellar metallicity translates to an [O/H] value of —1.13 dex 
(Mateo [19981 , or to a 12+log[0/H] value of 7.78, therefore 
is consistent with the metallicity-luminosity relation in Lee et 
al. ( 120071 1. Moreover, the stellar metallicity of SDIG is consistent 
with the stellar metallicity-luminosity relation in Cloet-Osselaer 
et al. (I20T21 . 

The derived properties of SDIG are listed in Table [2] Each 
row lists the following information: (1) the average SFR; (2) and 
(3) the current and recent SFRs, respectively (see below); (4) 
the average age of the stellar populations; (5) the astrated stellar 
mass; (6) the stellar mass for stars younger than 1 Gyr; (7) the 
current stellar mass, M* iCunent ; (8) the gas mass, computed as 
1.4xMm; (9) the baryonic mass, computed as M gas +M+ jCurren t; 
(10) the baryonic gas fraction, computed as M gas /Mb a ryonic; (11) 
the current V-band stellar mass-to-light ratio; (12) the current 
B-band stellar mass-to-light ratio; (13) the average metallicity 
based on the SFH results; (14) the stellar metallicity based on the 
RGB fiducial and Padova isochrones; (15) the stellar metallicity 
based on the RGB fiducial and Dartmouth isochrones. 
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Table 2. Derived properties of SDIG. 



Quantity 



Value 



<SFR> 

SFRcurreat 
J1 AV recent 
<T> 

M, 

M gas 

■Mgajj / Mbaryonic 

M 4 /L v 

M,/L B 

[M/H] SFH 

[M/H] Padova 

[M/H] 

Dartmouth 



1.3^ xl(T J M© yr l 
2.5* xl(T 3 Mq yr- [ 
2.7±0.5xl(T 3 Moyr 1 
6A + _\1 Gyr 

1 - 77 -B73 Xl ° 7 M © 
1.68 + ™ xl0 6 M o 

L50™\| xlO 7 M© 

3.5x"l0 7 M© 

5xl0 7 M 

0.7 

3.2M O /L 
1.4Mo/L 

-1.5 dex 

-1.5 dex 

— 1.6 dex 



We can compare the SFR in our youngest age bin to current 
SFR indicators like Ha and FUV emission. The current SFR is 
estimated from our data assuming it is equal to the average SFR 
in our first age bin (4-30 Myr), 2.5+°g xl(T 3 Mq yr -1 . Ha emis- 
sion is sensitive to SFR over the past 10 Myr but can be biased 
to lower values than the true current SFR when the SFR is less 
than w3xl0~ 3 M yr -1 (e.g., Lee et al. 120091 ) because it relies 
on ionising photons from O stars, which may not be fully sam- 
pled by the IMF in cases of low SFR. We detect of order a dozen 
stars bright enough and blue enough to be candidate OV stars, 
so stochastic effects are expected to be important. Despite this 
caveat, the current SFR we derive from the CMD modelling is 
consistent with the upper limit of 2xl0~ 3 Mq yr~ l that Heisler 
et al. (119971 1 derive based on their 3 cr detection of the Ha sur- 
face brightness across the continuum-subtracted image of SDIG. 
However, Bouchard et al. (2009) derive a lower current SFR of 
3.8xl0~ 5 Mq yr~' based on the Ha flux in two detected knots. 
Their SFR is just 3% of the value we find. This could be taken 
as evidence for a quiescent period of low SFR, or for random 
variance due to the low numbers of massive stars formed. 

We can estimate the number of O-type stars expected to be 
produced by the recent SFR using the stellar population synthe- 
sis code Starburst99 (Leitherer et al. 119991 ). As a rough estimate, 
we assume continuous star formation at a constant rate equal to 
the recent SFR (see Table 13, a Salpeter IMF, and solar metal- 
licity. The results provide us with the total number of ionising 
photons emitted during the most recent time step of 2 Myr, and 
over a longer period of 100 Myr. Using the ionising flux Qo for 
an 05 V star (Q = 10 49 22 s" 1 , M *37 M ; Martins, Schaerer & 
Hillier 2005 ), we estimate that of order 18 such stars would have 
been produced over 100 Myr given our CMD-derived SFR. Of 
these, only a small fraction would be expected to be observed 
today owing to their short lifetimes. As an additional way to es- 
timate the number of massive stars present, we employ a Monte 
Carlo approach to estimate the variance in massive star numbers 
given an overall low SFR. To that end, we assume a Salpeter 
IMF (Salpeter 119551) with lower and upper mass limits of 0.5 
to 150 Mq, and draw from this distribution until 2.7xl0 5 Mq 
total mass of stars are formed; this number corresponds to our 
recent SFR integrated over 10 8 yr. The process is repeated 100 
times, allowing us to calculate that on average 1216+40 stars 
with masses appropriate to O stars (here assumed to be 18- 
100 Mq) are formed over this time period. At any given moment 
we would expect to see ^120 O stars if the IMF was being fully 
sampled, assuming an average lifetime of 10 Myr. The fact that 



we observe about an order of magnitude less confirms that SDIG 
is in a "down" state of star formation lasting at least the duration 
of a typical O star lifetime. However, given the SFR determined 
over our youngest age bin (4-30 Myr) is higher than the life- 
time average for the galaxy and the high gas content of SDIG, 
it seems more likely that this represents a combination of short- 
term SFR fluctuation and stochastic effects due to IMF sampling 
than to any physically meaningful suppression of star formation 
on longer timescales. 

Because Ha is known to be problematic in cases of very low 
SFR, the FUV emission provides a better estimator of the recent 
SFR, sampling stars down to a few solar masses, and ages of 
order 100 Myr (Lee et al. 1201 II ). We use t he FUV asymptotic 
magnitude from Table 2 of Lee et al. 1201 II in conjunction with 
the FUV-SFR relation from Lee et al. ([2009; their §3) to estimate 
a recent SFR of 9.6xl0~ 4 Mq yr~' . For comparison, the SFR we 
derive from the CMD is 2.7±0.5 xlO 3 M yr -1 . We therefore 
find improved agreement between the recent SFR indicator and 
the CMD modelling, but a factor of 2-3 discrepancy remains, 
presumably due to the fact that a galaxy that is deficient in O 
stars will also be deficient in FUV flux. 



6. Summary and discussion 



SDIG is a gas-rich, low-luminosity dlrr (Cote et al. l 1997t Heisler 
et al. 119971 ). It is dominated by intermediate-age stars with an av- 
erage age t =6.4+ J'Jj Gyr, similar to the mean mass-weighted age 
of LG dlrrs (e.g., Orban et al. 120081 ). The average SFR is simi- 
lar to other Sculptor group dwarfs, listed in Weisz et al. ([201 lb; 
their Table 2), with ESO540-G030 and ESO540-G032 having 
the lowest and highest average SFRs, respectively. The non- 
detection of H II regions and low Ha-based SFR have been re- 
garded as indicative of SDIG experiencing a currently quiescent 
phase (Miller [19961 Heisler et al. [19971 Skillman et al. I2003at 
Bouchard et al. 2009). However, our SFH modelling shows ele- 
vated star formation for ages younger than 100 Myr, consistent 
with the star formation levels detected in the FUV emission. Our 
current data do not allow us to resolve star formation events last- 
ing ~10 8 Myr for ages older than ~1 Gyr. There are intriguing 
hints of a strong previous enhancement of SFR at ages of 600- 
1100 Myr, although the amplitude and timing of this event are 
model-dependent. Similar patterns of globally decreasing SFR 
during the past 1-2 Gyr have been observed in a number of LG 
galaxies with similar luminosity to SDIG, e.g. Leo A (Cole et 
al. l2007l ). 

SDIG is a member of the NGC 7793 subgroup located at the 
far-side of the Sculptor group (Jerjen et al. 1998. Karachentsev 
et al. 120061) . and NGC 7793 is identified as its main disturber, 
i.e. the neighbouring galaxy producing the maximum tidal ac- 
tion (Karachentsev et al. 2004). NGC 7793 is a spiral galaxy and 
stellar population analyses of several halo and disk fields detect a 
metallicity gradient (Vlajic et al. 1201 II ). as well as indications for 
substantial stellar radial migration (Radburn-Smith et al. 120121 ). 
The SFH of a halo field in NGC 7793 indicates that a burst of 
star formation has occurred within the last 20 Myr and with an 
amplitude three times that of the average SFR, most likely as- 
sociated with the spiral arm feature contaminating the halo field 
(Dalcanton et al. [20"T2l Figs. 1, 2, and 9 for NGC 7793). The de- 
projected distance between SDIG and NGC 7793 is 4 50 +34 jj kpc, 
estimated using the method in Karachentsev et al. J2004t . i.e. 
R 2 = D 2 + D^p - 2DDMDCOS0, where Dmd is the distance in 
Mpc of NGC 7793 adopted from Vlajic et al. J20TTT ). D is the 
distance in Mpc of SDIG, and is their angular separation, cal- 
culated using NASA/IPAC Extragalactic Database (NED). The 
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uncertainty in the deprojected distance takes into account the 
uncertainties in the distance of NGC 7793 and SDIG, as well as 
their angular projected separation as the minimum physical dis- 
tance between them. We note that many gas-rich irregular and 
transition dwarfs are similarly distant from the Milky Way, for 
example NGC 6822 or LeoT (McConnachie 2012) . Given the 
mass of NGC 7793 from Carignan & Puche (fl990T > the orbital 
timescale for SDIG around NGC 7793 would appear to be too 
long to be responsible for star formation activity on timescales 
shorter than a Hubble time, but other signs of disturbance are 
present. 

Cote et al. (2000) discuss how the low value of SDIG's 
maximum rotational velocity implies that random motions may 
significantly contribute to the dynamical support and may thus 
lead to the puffing-up of an initially thin disk. To that end, 
the twisted Hi disk of SDIG (Cote et al. 2000J may indicate 
an interaction having occurred in this dwarf's past. NGC 7793 
has been mapped with Chandra (Pannuti et al. 1201 11 1, and the 
results indicate an asymmetry in the detected X-ray sources. 
Pannuti et al. (1201 11 1 discuss an interaction between SDIG and 
NGC 7793 as a plausible way to explain this asymmetry in the 
X-ray sources, although they do not regard it as strong a pos- 
sibility. Radburn-Smith et al. (120121 ) discuss that the detected 
stellar radial migration in NGC 7793 may also be due to a past 
interaction event having occurred in NGC 7793. While the de- 
projected distance of SDIG shows its current location and its 
actual orbit remains unknown (e.g. Bellazzini et al. 1996]), the 
timescale of the star formation enhancements and the crossing 
timescale of SDIG around NGC 7793 suggest that the elevated 
SFRs in SDIG are not due to gravitational interactions between 
them. The SF enhancements and disturbed H i kinematics seen 
in SDIG may be due to the effect of gravitational interactions 
with a lower luminosity, gas-poor dwarf galaxy companion that 
has thus far remained elusive, as such tidal interactions between 
dwarf galaxies now start to become revealed (e.g., Martinez- 
Delgado et al. 120121 Rich et al. 120121 Hunter et al. [19981 . Tidal 
interactions between dwarf galaxies may also provide an impor- 
tant link to understanding the evolution and characteristics of 
dwarf galaxies in the isolated far-outskirts of the LG, as the case 
of the transition-type dwarf galaxy VV124 shows (e.g., Kirby et 
al. 2012). 

It is interesting in this respect that based on photometric ar- 
guments alone, Heisler et al. (1997) conclude that it is not pos- 
sible to rule out that SDIG-like systems may evolve to dSphs. 
It is thus intriguing that the off-center local sites of star forma- 
tion observed in the young MS stars of SDIG, as well as other 
LG d warfs ( Dohm-Palm er et a l. [19971 GaUagher e t al. [19981 
Mateo [19981 Cole et al. IT9991 McConnachie et al. 120061 are 
observed in the majority of the transition-type dwarfs in the 
Sculptor group (Lianou et al. 120121 Jerjen & Rejkuba 2001). 
Environmental effects on gas-rich dwarfs provide a possible av- 
enue to their evolution into gas-poor systems (e.g., Mayer 2010; 
and references therein), and such environmental effects seem 
to affect their recent star formation activity (e.g., Bouchard et 
al. 120091 Thomas et al. r20"T0l Weisz et al. 1201 lbb . However, the 
detection of a recent modest upturn in SFR, and hints of a simi- 
lar episode ^0.6-1 Gyr ago is further indication that small galax- 
ies are able to sustain complex star formation histories without 
strong tidal interactions from larger neighbours. 
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